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Abstract 

We recall the non-Abelian Stokes theorem for the Wilson loop in the Yang -Mills 
theory and discuss its meaning. Then we move to 'gravitational Wilson loops', i.e. to 
holonomies in curved d = 2,3,4 spaces and derive 'non-Abelian Stokes theorems' for 
these quantities as well, which are similar to our formula in the Yang-Mills theory. In 
particular we derive an elegant formula for the holonomy in the case of a constant- 
curvature background in three dimensions and a formula for small-area loops in any 
number of dimensions. 



1 Introduction 

One of the main objects in the Yang-Mills theory as well as in gravity is the parallel trans- 
porter along closed contours, or holonomy. In Yang-Mills theory it is conventionally called 
the Wilson loop; it can be written as a path-ordered exponent, 

I r Ht^ 1 

W r = — TrP exp i 6 dr — A a T a , (1) 
d[r) J dr p 

where x M (r) with < r < 1 parametrizes the closed contour, A ^ is the Yang-Mills field (or 
connection) and T a are the generators of the gauge group in a given representation r whose 
dimension is d(r). In curved Riemannian spaces the 'gravitational Wilson loop' or holonomy 



for d- dimensional vectors can be also written as a trace of the path-ordered exponent of the 
connection, this time of the Christoffel symbol, 



One can also consider parallel transporters of spinors in curved background: in this case the 
holonomy is defined not by the Christoffel symbols but by the spin connection which is not 
uniquely determined by the metric tensor, see the precise definitions below. 

The Yang-Mills Wilson loop is invariant under gauge transformations of the background 
field A^; the gravitational Wilson loop is invariant under general coordinate transformations 
or diffeomorphisms, provided one transforms the contour as well. 

It is generally believed that in three and four dimensions the average of the Wilson loop in 
a pure Yang-Mills quantum theory exhibits an area behaviour for large and simple contours 
(like flat rectangular). This should be true not for all representations but those with W-ality' 
nonequal zero; in the simplest case of the SU (2) gauge group these are representations with 
half-integer spin J. 

One of the difficulties in proving the asymptotic area law for the Wilson loop in half- 
integer representations (and proving that in integer representations it is absent) is that the 
Wilson loop is a complicated object by itself: it is impossible to calculate it analytically 
in a general non-Abelian background field. Meanwhile, it is sometimes easier to average a 
quantity over an ensemble than to calculate it for a specific representative. However, in case 
of the Wilson loop the path-ordering is a serious obstacle on that way. 

A decade ago we have suggested a formula for the Wilson loop in a given background 
belonging to any gauge group and any representation [|TJ. In this formula the path order- 
ing along the loop is removed, but at the price of an additional integration over all gauge 
transformations of the given non-Abelian background field, or, more precisely, over a coset 
depending on the particular representation in which the Wilson loop is considered. Fur- 
thermore, the Wilson loop can be presented in a form of a surface integral 0, see the next 
section. We call this representation the non-Abelian Stokes theorem. It is quite different 
from previous interesting statements [B[ f|, [], also called by their authors 'non-Abelian 
Stokes theorem' but which involve surface ordering. Our formula has no surface ordering. A 
classification of 'non-Abelian Stokes theorems' for arbitrary groups and their representations 
has been given recently by Kondo et al. J7| who used the naturally arising techniques of flag 
manifolds. 

Though these formulae usually do not facilitate finding Wilson loops in particular back- 
grounds, they can be used to average Wilson loops over ensembles of Yang-Mills configura- 
tions or over different metrics, and in more general settings, see. e.g. || [5], |7|, |H| . 

The main aim of this paper is to present new formulae for the gravitational holonomies 
in curved d = 2,3,4 spaces: they are similar to our non-Abelian Stokes theorem for the 
Yang-Mills case. We get rid of the path ordering in eq. (|2|) and write down the holonomies 
as exponents of surface integrals. Instead of path ordering we have to integrate over certain 
covariantly unit vectors (in case of d = 3) or covariantly unit (anti) self-dual tensors (in case 
of d = 4). Remarkably, these formulae put parallel transporters of different spins on the 
same footing. In particular, holonomies for half-integer spins are presented in terms of the 
metric tensor (and its derivatives) only but not in terms of the vielbein or spin connection. 
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Apart from purely theoretical interest we have a practical motivation in mind. Quite 
recently we have shown, both in the continuum and on the lattice, that the SU(2) Yang- 
Mills partition function in d — 3 can be exactly rewritten in terms of local gauge-invariant 
quantities being the six components of the metric tensor of the dual space [11]. This rewriting 
may be useful to investigate the spectrum and the correlation functions of the theory directly 
in a gauge-invariant way, but it is insufficient to study the interactions of external sources 
since they couple to the Yang-Mills potential and not to gauge-invariant quantities. The 
present paper demonstrates, however, that a typical source, i.e. the Yang-Mills Wilson loop 
can be expressed not only through the potential (or connection) but also through the metric 
tensor which is gauge-invariant. Thus, not only the partition function but also Wilson loops 
in the d = 3 Yang-Mills theory can be expressed through local gauge-invariant quantities. 
We leave a detailed formulation of the resulting theory for a forthcoming publication. 

Though the main content of the paper are the non-Abelian Stokes theorems for holonomies 
in 3 and 4 dimensions, we have added three short sections with relevant material. We add 
for completeness the Stokes theorem in two dimensions, compute the holonomy in a special 
case of constant curvature with cylinder topology in three dimensions and give a general 
formula for the 'gravitational Wilson loop' for small loops in any number of dimensions. 



2 Non-Abelian Stokes theorem in Yang— Mills theory 

Let t parametrize the loop defined by the trajectory x M (r) and A(r) be the tangent compo- 
nent of the Yang-Mills field along the loop in the fundamental representation of the gauge 
group, A(t) = A a ^t a dx^ /dr, Tr(t a i 6 ) = \5 ab . The gauge transformation of A{r) is 

A(r) - S(t)A(t)S-\t) + iS{r)^S-\r). (3) 

Let Hi be the generators from the Cartan subalgebra (i = 1, r; r is the rank of the gauge 
group) and the r-imensional vector m be the highest weight of the representation r in which 
the Wilson loop is considered. The formula for the Wilson loop derived in ref. is a path 
integral over all gauge transformations £(t) which should be periodic along the contour: 

W r = J DS(t) exp i J dr Tr [m t Hi (SAS' 1 + iSS' 1 )] . (4) 

Let us stress that eq. (f|) is manifestly gauge invariant, as is the Wilson loop itself. For 
example, in the simple case of the £77(2) group eq. (|) reads: 

Wj = J DS(t) exp i J J dr Tr [t 3 (SA& + z££ f )] (5) 

where t% is the third Pauli matrix and J — |, 1, |, ... is the 'spin' of the representation of 
the Wilson loop considered. 

The path integrals over all gauge rotations (|],[5J) are not of the Feynman type: they do 
not contain terms quadratic in the derivatives in r. Therefore, a certain regularization is 
understood in these equations, ensuring that £(r) is sufficiently smooth. For example, one 
can introduce quadratic terms in the angular velocities iSS^ with small coefficients eventually 
put to zero; see ref. |lj] for details. In ref. fl]] eq. (§]) has been derived in two independent ways: 
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i) by direct discretization and ii) by using the standard Feynman representation of path 
integrals as a sum over all intermediate states, in this case that of an axial top supplemented 
by a Wess-Zumino' type of the action. Another discretization but leading to the same 
result has been used recently by Kondo 0. A similar formula has been used by Alekseev, 
Faddeev and Shatashvili ]12| who derived a formula for group characters to which the Wilson 



loop is reduced in case of a constant A field actually considered in [TJ]. In ref.p3[ eq. (f|) 
has been rederived in an independent way specifically for the fundamental representation 
of the SU(N) gauge group. Finally, another derivation of a variant of eq. ([5]) using lattice 
regularization has been presented recently in ref. Jl4j] . 

The second term in the exponent of eqs. is in fact a Wess-Zumino '-type action, 

and it can be rewritten not as a line but as a surface integral inside a closed contour. Let 
us consider for simplicity the SU(2) gauge group and parametrize the SU(2) matrix S from 
eq. fl5|) by Euler's angles, 

S = exp(ryr 3 /2) exp(z/5r 2 /2) exp(iar 3 /2) 

(6) 



cos I e 2 sin | e 2 
— sin | e 2 cos ^ e 2 



The derivation of eq. (||) implies that S(t) is a periodic matrix. It means that a ± 7 and (3 
are periodic functions of r, modulo Air. 

The second term in the exponent of eq. (||) which we denote by $ is then 

$ = J dr Tr(r 3 iS&) = J dr d(cos (3 + 7) 

= Jdr[a(cosp- 1) + (d + 7)] = JdTa(cosf3- 1). (7) 

The last term is a full derivative and can be actually dropped because a + 7 is 47r-periodic, 
therefore even for half- integer representations J it does not contribute to eq. Notice 
that a can be 27r-periodic if 7 (which drops from eq. ([?[)) is 2tt, 6tt, . . .-periodic. If a(l) = 
a(0) + 2irk, a{r) makes k windings. Integration over all possible a(r) implied in eq. @ can 
be divided into distinct sectors with different winding numbers k. 
Introducing a unit 3-vector 

n a = -Tr (S , r a 5 , V 3 ) = (sin (3 cos a, sin j3 sin a, cos/5) (8) 

we can rewrite $ as 

$ = - / drda e abc e ij n a d i n b d j n c , i,j = r,a, (9) 



2 

where one integrates over any surface spanned on the contour (we shall call it a 'disk'), and 
n or a and (3 are continued inside the disk without singularities. We denote the second 
coordinate by a such that a = 1 corresponds to the edge of the disk coinciding with the 



contour and a = corresponds to the center of the disk. See ref. |14[ for details on the 
continuation inside the disk. 



4 



Let us note that if the surface is closed or infinite the r.h.s. of eq. 
topological charge of the n field on the surface: 



Q = —j dadre abc e ij n a d i n b d j n c 



is the integer 



(10) 



Eq. (|9|) can be also rewritten in a form which is invariant under the reparametrizations 
of the surface. Introducing the invariant element of a surface, 



d 2 S fMU = dadr 
one can rewrite eq. @ as 

$ = 

We get for the Wilson loop [ffl: 



' dx^ dx u dx u dx^ " 



(IV 



dr da dr da 



d 2 S^e abc n a d„n b d„ 



d(Area), 



■nr. 



Wj = j Dn(a, r) exp iJ J dr(A a n a ) + y / d 2 S"" 'e abc n a d^n b d } 



(12) 



(13) 



The interpretation of this formula is obvious: the unit vector n plays the role of the 
instant direction of the colour 'spin' in colour space; however, multiplying its length by J 
does not yet guarantee that we deal with a true quantum state from a representation labelled 
by J - that is achieved only by introducing the 'Wess-Zumino' term in eq. QIB"!): it fixes the 
representation to which the probe quark of the Wilson loop belongs to be exactly J. 

Finally, we can rewrite the exponent in eq. flT3"| ) so that both terms appear to be surface 
integrals 0: 



W 



jDn(a,r) exp y J d 2 S^ (-J£, 



n a + e abc n a ^ji^jc 



(14) 



where _D" fc 

h 1 



d,,5 ab + e ach A c „ is the covariant derivative and F, a , 



A a — f) Aa \ abc Ab Ac 

is the field strength. Indeed, expanding the exponent of eq. (|1~4| ) in powers of one observes 
that the quadratic term cancels out while the linear one is a full derivative reproducing the 
A a n a term in eq. (|T3"D; the zero-order term is the 'Wess-Zumino' term (^) or ([7j). Note that 
both terms in eq. ([U]) are explicitly gauge invariant. We call eq. (0) the non-abelian Stokes 
theorem. We stress that it is different from previously suggested Stokes-like representations 
of the Wilson loop, based on ordering of elementary surfaces inside the loop 0, |], [5|, || . For 
a further discussion of eq. QT3D see [TJJ . 

Let us briefly discuss gauge groups higher than SU (2): for that purpose we have to return 
to our eq. (|). Eq. @) is valid for any group and any representation. It is easy to present it 
also in a surface form, see ref. ||14|| . Actually, eq. (|J) depends not on all parameters of the 
gauge transformation but only on those which do not commute with the Cartan combination 
TC r = rriiHi. In the SU{2) case one has ni^Hi = Jr 3 , J = 1/2, 1, 3/2, . . ., since £77(2) is of 
rank 1, and there is only one Cartan generator. Therefore, in the £77(2) case one integrates 
over the coset SU(2)/U(1) for any representation; this coset can be parametrized by the n 
field as described above. 



5 



In case of higher groups the particular coset depends on the representation of the Wilson 
loop. For example, in case the Wilson loop is considered in the fundamental representation 
of the SU(N) group the combination mj-ffj is proportional to one particular generator of the 
Cartan subalgebra, which commutes with the SU (N — 1) x U(l) subgroup. [In case of SU (3) 
this generator is the Gell-Mann Ag matrix or a permutation of its elements.] Therefore, the 
appropriate coset for the fundamental representation of the SU (N) group is SU(N) /SU (N— 
1) /U{1) = CP 1 ^^ 1 . A possible parametrization of this coset is given by a complex iV-vector 
u a of unit length, u^u 01 = 1. To be concrete, the Cartan combination in the fundamental 
representation can be always set to be mjifj = diag(l, 0, . . . , 0) by rotating the axes and 
subtracting the unit matrix. In such a basis u a is just the first column of the unitary matrix 
S while li* is the first row of S^ . Unitarity of S implies that u^ a u a = 1. 

In this parametrization eq. (H) can be written as 



W, 



SU{N) 



J Du Du^ 5(ulu a - 1) exp i J dr ^£u{ (zV^ vP, (V„) J = d^- iA^ {t% . 

(15) 

Using the identity, 

(V;w) f (VjU) +M t V i V i 



fund 



(16) 



we can present eq. fllBD in a surface form: 



J DuDv)5{\u\ 2 - 1) expiJdS^ ]^{u ] F^u) + i(V » f {S7 v u 



;i7) 



where F^ is the field strength in the fundamental representation. Eq. ([17]) has been first pub- 
lished in ref.|13| however with an unexpected overall coefficient 2 in the exponent. Eq. (IT/ 



presents the non-Abelian Stokes theorem for the Wilson loop in the fundamental represen- 
tation of SU(N). In the particular case of the SU(2) group transition to eq. (|H]) is achieved 
by identifying the unit 3-vector: n a 



u' a (T a )puP where 



u 



cos ^ e 2 
sin | e 1 2 



2i v)d T u 

It should be mentioned that the quantity 



d(cos P — 1) + (a + 7). 



J dadreij idiv) a djU a = 2nQ 



(19) 

_1 model. For closed 



appearing in eq. (|17D is the topological charge of the 2-dimensional CP N ~ 
or infinite surfaces Q is an integer. 

In case the Wilson loop is taken in the adjoint represention of the SU (N) gauge group the 
combination m^ii/j is the highest root. Only group elements of the form exp(zajifj) commute 
with this combination, belonging to the maximal torus subgroup ^/(l)^ -1 . Hence, in case of 
the adjoint representation one in fact integrates over the maximal coset SU(N)/ [/(l)^ -1 = 

F N- 



i.e. over flag variables |T5 
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3 'Gravitational Wilson loops' 



An object similar to the Wilson loop of the Yang-Mills theory exists also in gravity theory. 
It is the parallel transporter of a vector on a Riemannian manifold along a closed contour, 
else called a holonomy. The holonomy is trivial if the space is flat but becomes a non-trivial 
functional of the curvature in case it is nonzero. In the remaining sections we shall present 
new formulae for the parallel transporters on d = 2, 3, 4 Riemannian manifolds. 



Let us first remind notations from differential geometry. We use |L6[ as a general reference 
book. 

Let g^ v = g Vil (//, v — 1, . . . d) be the covariant metric tensor, with the contravariant g^ v 
being its inverse, g fl u9 UK = The determinant of the covariant metric tensor is denoted by 
g. The Christoffel symbol is defined as 



;<A 



T» K = g» x T XyUK = 9 —{d v g XK + d K g Xv - d x g v 



2g ' 



The action of the covariant derivative on the contravariant vector is defined as 



(v p )X = («Va + r;> A . 



The commutator of two covariant derivatives determine the Riemann tensor: 



(20) 



(21) 



[V P V CT ]^ - R\ pu - g KR R K ' Xpa - d p T K aX - d ff Y K pX + T K pT T T aX - r^ r r^ A . 
A contraction of the Riemann tensor gives the symmetric Ricci tensor, 



(22) 



R 



Act 



R\ 



Akct' 



Its full contraction is the scalar curvature: 



r; 



R — R\a9 Xcr — R1- 



DK _Acr 
K \pa9 ■ 



(23) 



(24) 



The parallel transporter of a contravariant vector along a curve x M (r) is determined by 
solving the equation, 

dx' 1 



dr 



(V M )> A (r) = 0, 



whose solution can be written with the help of the evolution operator, 



W g (t 



v x (0) 



(25) 



(26) 



where t> A (0) is the vector at the starting point of the contour and v x (t) is the parallel- 
transported vector at the point labelled by r. The evolution operator can be symbolically 
written as a path-ordered exponent of the Christoffel symbol: 



W 



Gi 



P exp 



dT——T p 
dr 



(27) 



J A 



We define the 'gravitational Wilson loop' as the trace of the parallel transporting evolution 
operator along the closed curve x^(t) with x^(l) = x^(0): 
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rG 



1 r 



' vector j 



(2J 



This quantity is diffeomorphism- invariant: if one changes the coordinates x M — > x'^(x) 
the metric is transformed, but if one changes the contour accordingly, that is x M (r) — > 
x' M (a;(r)) the gravitational Wilson loop or the holonomy remains the same. In this respect 
the gravitational holonomy is different from the Yang-Mills Wilson loop which is invariant 
under gauge transformation, without changing the contour. 

The parallel transporter of a covariant vector is given by the transposed matrix; its trace 
coincides with that of the contravariant vector. 

4 Relation of gravity quantities to those of the Yang- 
Mills theory 

We shall show now that the 'gravitational Wilson loop' is not just analogous but directly 
expressible through the Yang-Mills Wilson loops of the SU(2) group. To that end we 
introduce the standard vielbein e A and its inverse e Afi such that 

4<£ = 9i*> e A e B » = 5 AB , e A ^ = g^, det e A = y/g. (29) 
Let us decompose the vector experiencing the parallel transport in vielbeins, 

v x = c A e AX , the reciprocal being c A = e A v K , (30) 
and put it into eq. (|25| ) defining the parallel transport. We have 



dx^ ,_ - k a ax dx^ 



= ^_(V M )^c A e AA = ^- e AK d^c A + c A (d^e AK + T K x e AX ) = ^-e BK (d^5 BA + uj ba )c a , 



dr v dr 



where we have introduced the spin connection, 



dx^ 



dr 



(31) 



<*f = ~< A = \e AK {d,e B K - d K e B ) - \e B \d»e A - d K e A ) - \e A «e BX e<;{d K e c x - d x e% (32) 

and used the fundamental relation: 

d,e AK + r; x e AX = -cu AB e B ", (33) 
d,e A - I>f = -ufe B . (34) 
One can introduce the SO(d) 'field strength', 

K B = % + «>» 9 » + = d ^ AB ~ + - <"W* (35) 

related to the Riemann tensor as 
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The above material is common for any number of dimensions. To proceed further we 
need to consider separately the cases d = 3 and d — 4. The case d — 2 is considered in 
section 6. 



4.1 (1=3 

In three dimensions one can immediately identify the spin connection with a Yang-Mills 
field in su(2), 



abc .ab 
— c lV: . 

2 



(37) 



Speaking about three dimensions we shall denote Lorentz indices by i,j, ... = 1, 2, 3 and the 
flat triade indices by a, b, ... = 1, 2, 3. Recalling the generators in the J = 1 representation, 



Crridl 



we can rewrite the last parenthesis in eq. (|31~D as 



l€ 



'dfrpf 



(38) 



\C t rr\c\ab 



+ uf = di5 ab - %A\{T 



(A 



i ab 



(39) 



which is the standard Yang-Mills covariant derivative in the adjoint representation. In the 
fundamental (spinor) representation the Yang-Mills covariant derivative is 



(VOS = dM -iAU-\ = M + -cj. 



a' 



.ab 



a a 



<x,P = l, 2, 



(40) 



which coincides with the known expression for the covariant derivative in the spinor repre- 
sentation in curved space. 

The standard Yang-Mills field strength is directly related to that of eq. Q35|) : 



F% = diA) - djAI + e ahc A\A) 



- abc rbc 



(41) 



Hence from eq. (|36|) one has 



e abc FZelef = R, 



ij^k^i — J-Ujki- (42) 

Let us consider the parallel transporter of a 3-vector in curved space, as defined by 
eq. (|25|) . According to eqs. ( |31~1 , |3~9D solving eq. is equivalent to solving the Yang-Mills 
equation for the parallel transporter, 



dx l 



(A) ab c b 



0. 



(43) 



whose solution is 



c a (r) 



YM i 



ab 



c b (0), 



YM, 



T 



ab 



P exp i Idr A C ;T L 

' dr 



ab 



(44) 



where the subscript "1" refers to that the path-ordered exponent is taken in the J 
representation. The parallel transport of a contravariant vector is therefore 
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(r) 



c a (T)e ak (T) 



e ak (r) 



YM, 



T 



ab 



ef (0)^(0), 



(45) 



from where we immediately get the needed relation between the 'gravitational' and Yang- 
Mills parallel transporters: 



G, 



T 



ab 



bl 



(0). 



(46) 



The relation becomes especially neat for the Wilson loops, i.e. for the traces of parallel 
transporters along closed contours. Since for the closed contour the vielbeins at the end 
points are the same, e%{l) = e£(0), we get 



G 

vector 



G 



W{ 



YM 



YM 



(47) 



In a similar way one can show that the same equation holds true for the gravitational parallel 
transporter of covariant vectors and, more generally, for parallel transporters of any integer 
spin J. In this case the Yang-Mills Wilson loop should be taken in the same representation 
as the gravitational one: 



W 



G 



w 



YM 



J -rrj ■ (48) 

In eq. fl48|) it is understood that the r.h.s. is expressed through the Yang-Mills field 
equal to the spin connection according to eq. (ff7|), while the l.h.s. is expressed through 
the Christoffel symbols, that is through the metric. It should be stressed that the spin 
connection is defined via the vielbein which is not uniquely determined by the metric tensor. 
Nevertheless, the Wilson loop, being a gauge-invariant quantity, is uniquely determined by 
the metric tensor and its derivatives. This is the meaning of eq. (48 



For half-integer J there is no way to define the parallel transporter other than through 
the spin connection. Nevertheless, as we show in section 8 where we present the holonomy 
for any spin in a surface form, the 'gravitational Wilson loop' is also expressible through the 
metric tensor and its derivatives, even for half-integer spins. 



4.2 d=4 

In four Euclidean dimensions the rotational group is £0(4) = SU{2) x 577(2), therefore all 
irreducible representations of 5*0(4) can be classified as (Ji, J 2) where Ji j2 = 0, ~, 1, ... are 
the representations of the two 577(2) subgroups. For example, the 4-vector representation 
whose parallel transporter has been considered in the beginning of this section, transforms 
as the (~, |) representation of SU{2) x 577(2). 

Because of this, it is convenient to decompose the spin connection u AB into self-dual and 
anti-self-dual parts using 't Hooft's 77 and fj symbols. They are defined as 

if AB = ^Tra a (a A+ a B - - a B+ a A ~), a A± = (±z<7,l), (49) 

f] aAB = ^Tia a {a A -a B + -a B -a A+ ). (50) 



We use capital Latin characters to denote flat 4-dimensional vierbein indices, A, B 



1, 2, 3, 4, while a,b, . . . = 1,2, 3; a a are the three Pauli matrices. The spin connection ui AB 
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transforms as a 6- dimensional representation of 5*0(4), which can be decomposed as a sum 
(1, 0) + (0, 1) of adjoint representations of the two SU(2) subgroups. We write: 



CO 



AB 



-a „aAB a -aAB 



P^V 



The 50 (4) 'field strength' (|35|) is then decomposed accordingly: 



AS 



-aAB 



(51) 



(52) 



where 



KM 



fj, 1 - " fi ' v 
/./'/ ' ~fJ.Pi; ^uPfi I'lil'r 



(53) 
(54) 



are the usual Yang-Mills field strengths of the SU(2) Yang-Mills potentials ir® and p a . It 
should be stressed that 6 ■ 4 = 24 variables ui AB equivalent to 2 ■ 3 ■ 4 = 24 variables it®, are 
defined by only 4 • 4 = 16 tetrades e A via eq. (|32"D, so that not all of them are independent. 
Contracting eq. ([36]) with the 7],fj symbols we get: 



K(p) 



2 
1 

2 r 



V aAB e AK e BX R K 



-aAB Ak BAn 



(55) 
(56) 



Let us now return to the parallel transporter of a 4-vector. As shown in the beginning 
of this section, finding it is equivalent to solving the equation 



dr 



o %AB | AB\ „B 



0. 



Let us present the 4-vector c A as a combination of two spinors, 
Putting it into eq. (|57| ) and decomposing u AB as in eq. ( |51~D we get: 



9„ 



5 ^ " 



aAB 



+ p a u n 



-aAB 



1,2. 



0. 



(57) 



{hi 



(59) 



Using the definition of the r/-symbols (|49| , |50D it is easy to check that this equation is satisfied 
provided the spinors x, satisfy 



dx^ 
dr 
dx^ 
dr 









(?: 






d s 


-*< 


)1 










a" 




ra: 

Of) 


-i Pi 


'J 



or 



0. 



dx» 
dr 



0. 



(60) 
(61) 
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Expressions in square brackets are identical to the Yang-Mills covariant derivatives, the role 
of the Yang-Mills potentials played by and p a , respectively. Eqs.(]60|, |61]) define the Yang- 
Mills parallel transporters in the fundamental representation. Their solution can be written 
as evolution operators, 



[W^rXx^O) or xl(r) = X \(0) [W^ (r) 
[W p (r)f s ^(0), 







fdr 


dx^ 




p 


exp i 


dr *» 


Y 






[dr 


dxf 


a a ~ 


p 


exp i 


dr P » 


T 



(62) 
(63) 

(64) 
(65) 



J i 



Returning to the 4-vector c A ([58]) we see that its evolution is determined by 

= [WW(r)] AB c B (0), 



c A (r) 



[Wvector(r)] 



AB 



Tr 



W^{r)a A+ W p {r)a 



B- 



(66) 



Let us take a closed contour and take the trace of the evolution operator. The 'gravitational 
Wilson loop' for a 4-vector is then 



WFi is 



\e A %l)[W vectoI (l)] AB e B K 

-Tr W* ■ -Tr W. 
2 2 



\ [w vcctoI (i)] AA 



(67) 



Its generalization to the holonomy in arbitrary representation (Ji, J 2) is obvious: 



wi-w p j2 , 



2TTT Tr ^^ P - 



(6t 



Thus, in curved d = 4 space the holonomy in the ( Ji, J 2 ) representation is equal to a product 
of two Yang-Mills Wilson loops where the role of Yang-Mills potentials is played by self-dual 
(7T°i) and anti-self-dual parts of the spin connection. In section 9 we show that both W n 
and W p can be written in terms of the metric tensor. 



5 Small Wilson loops 

For small-area contours the 'gravitational Wilson loop' can be expanded in powers of the 
area. The most straightforward way to do it is to use the path-ordered form of W G as 
given by eq. (p7|). We take a square contour of size a x a lying in the 12 plane, and expand 
the path-ordered exponent in powers of a. After some simple algebra we obtain the first 
nontrivial term of that expansion which happens to be 0(a 4 ): 



G 



d 



G 

vector 



1 + J R * 



A12 rL k12 



1 - 



2(Asy u (Asy 

4d 



-R^uR P ^,'9 Kp 9 Xa , (69) 
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where (AS')'"' is the element of the surface lying in the (/iz/) plane. Notice that the first 
correction to the holonomy is negative-definite. It should be stressed that the first-order 
term in AS is, generally, present in the expansion of the parallel transporter, however it 
vanishes when one takes the trace owing to the identity R K Kfll/ = 0, so that the expansion for 
the trace starts from the (AS 1 ) 2 term. 

In d = 3 eq. (|69"D can be further simplified because the Riemann tensor is directly related 
to the Ricci tensor: 

R 

Rijkl = Rik9jl — RiWjk + RjWik — Rjkgil + -^(.9iWjk — gik9jl)- (70) 

Since Riemann tensor is antisymmetric inside both pairs of subscripts we can replace 

g km g ln - \{g km 9 ln - g kn g lm ) = ^e m r^ 9ij . (71) 

We introduce the dual element of the surface, 

AS pq = e pqr AS r , (72) 

and have 

e kU e^ r R klpq = -4 (V - ^RgA , (73) 

which as a matter of fact is the Einstein tensor. For the parallel transporter of arbitrary 
spin J the factor "2" in the numerator of (^) should be replaced by J (J + 1). 
Combining all the factors we obtain 

Wf = l- 2 J( ^ + 1} (iT - l -Rg ir ^j 9lJ (V - \Rg^ AS r AS S . (74) 

This is our final expression for the trace of the parallel transporter of spin J in a curved 
d = 3 space, for small loops. Notice that eq. ([74]) is invariant under diffeomorphisms. 



6 Gravitational Wilson loop in two dimensions 

In curved d = 2 space the trace of the parallel transporter along a closed loop can be 
computed exactly for any metric and presented in the form of a 'Stokes theorem'. The result 
is related to the Gauss-Bonnet theorem and is generally known: we present it here for the 
sake of completeness. 

The key observation is that in two dimensions the spin connection (|32| ) has only one 
component, 

ujf = e ab Ui. (75) 

In this section all indices assume only two values 1,2. According to eq. (|31|) the parallel 
transporter of a vector is determined by the equation 

_ ^ e ab c b = (76) 
dr dr 
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whose solution is 



c°(r) = W« b (r)c b (0), W ab (r) = ( COS ^ T f \ ^\ ) , 7 (r) = f ' dr^. 
K J v ; v V - sm 7(7) cos7(r) ^ ' !K J Jo dr 

(77) 

According to the general theorem of section 6 the gravitational Wilson loop is equal to the 
Yang-Mills one, and we obtain 

W? = ^W aa {l) = cos$ (78) 

where 

$ = 7 (1) = f Q dr^- Ui = y dx l e ab uf. (79) 

This formula is not fully satisfactory yet since the holonomy is expressed through the spin 
connection and not through the metric. It can be achieved if we apply the Stokes theorem 
and write eq. (|7"9|) in a surface form. We have 

$ = lj 'dSe^diuf (80) 

where dS is the element of the surface spanned on the contour. Introducing the field strength 
related to the Riemann tensor, 

F$ = d t uf - dpf + ufuf - ufuf = R% e a k e b , e ab e a k e b = e kl y/g, (81) 
and noticing that in d = 2 the commutator term is zero, we rewrite eq. flSOf ) as 

$ = \jdS^R, Wf = cos$, (82) 

where R = (1/2) e lj e^i R^j is the scalar curvature. It is gratifying that the holonomy is ex- 
pressed through the Einstein-Hilbert action, known to be a full derivative in two dimensions. 
Needless to explain, eq. (|82|) is diffeomorphism- invariant. 

In d = 2 there is essentially only one component of the Riemann tensor, 

#1212 = -Rg (83) 

(see ||16||). Taking it into account it is easy to check that for small areas the expansion of 
eq. ( |32"D gives the same result as eq. fl69|) written for small loops. 

7 An example of big loops: constant curvature back- 
ground in d = 3 

In three dimensions the Riemann tensor is expressible through the Ricci tensor, see eq. ([70]). 
Because of it, diffeomorphism-invariant information about curved spaces is fully contained 
in the three eigenvalues of the symmetric Ricci tensor, 



14 



i2$ = A<5$, (84) 

the scalar curvature being the sum of the three, R — \\ + A 2 + A3. For example, the 
de Sitter S 3 space corresponds to Ai = A2 = A3 = R/3 = const. In this section we would like 
to consider another special case of constant curvature, namely the cylinder space S 2 x R, 
characterized by Ai = A2 = R/2 = const, A3 = 0. We shall see that the parallel transporter 
in such spaces can be computed for any form of the contour and any metric and that the 
gravitational Wilson loop is given by an elegant formula. 

A general metric can be considered as induced by 6 external coordinates w (xi, X2, #3): 

9ij = d iW A d jW A , A = l,...,6. (85) 

In the special case of the cylinder space S 2 x R it is sufficient to use only four external 
coordinates w a (a = 1, 2, 3) and w 4 , subject to the constraint 

EK) 2 = |. (86) 

0=1 K 

An example of such external coordinates is given by 

F~2 T l,2,3 [W 

w 1A3 (x) = J- - , w\x) = J- lnr, (87) 



leading to the metric tensor 



3 



2 1 / 2 \ 2 1 

9 - = R^ 5 ^ ^={r) ^ (88) 

A simple calculation using formulae from section 3 shows that this metric indeed gives one 
eigenvalue of the Ricci tensor zero and the other two equal to a constant R/2. Since the 
eigenvalues of the Ricci tensor are diffeomorphism- invariant, a general change of coordinates, 
x % _> y % (x), in eq. (jB?|) result in the same eigenvalues. Therefore, the most general description 
of the cylinder spaces S 2 x R is given by 

„, \ [2y a (x) . [J 1 . . .. 2 diy a diy a . , 

1 — [2y 1 ijk diy a djy b d k y c 1 [R ijk a b 4 . . 

y/g = *e abc ^ = -e t3K e abc d i w a d j w°w c d k w*, (90) 

where y a (x) are three arbitrary functions of coordinates x i . Notice that gij is given by a 
product of two matrices Mf = diy a /\y\ and hence ^fg is itself a determinant (of the matrix 
M). 

Our aim is to calculate the Wilson loop for any contour in any metric (|89|) corresponding 
to the cylinder spaces. We shall make use of the fact that the Wilson loop is diffeomorphism- 
invariant. If we compute it for a general contour in some metric representing cylinder spaces, 
the most general case is recovered by diffeomorphisms of both the contour and the metric. 
We shall start with the specific metric given by eqs. flS7] , |8~8|) . 
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Given the metric tensor fl5?f), we construct a vielbein corresponding to it. This is, of 
course, not unique but any choice of the vielbein will suit us. We choose 




e: = l7;C ''I'" 9iy (91) 



Given the vielbein we construct the spin connection from its definition ( p2|) or the Yang-Mills 
field and obtain 

A a % = -^e ab Vf c = e aij 4, (92) 

which happens to be the field of the Wu-Yang monopole; the scalar curvature R has dropped 
from the spin connection. According to the theorem of section 4 the gravitational Wilson 
loop is equal to the Yang-Mills Wilson loop, provided the Yang-Mills potential A? is the 
spin connection of the metric under consideration. Therefore, all we have to do is to compute 
the Wilson loop in the field of the Wu-Yang monopole, but for a general contour. 

This task is easily solvable if we make use of another invariance, this time it is gauge 
invariance of the Wilson loop. It is well known that the Wu-Yang monopole in the hedgehog 
gauge (|92D can be transformed to the string gauge where the potential has only one nonzero 
component along the third colour axis (plus a Dirac string). In this gauge the Yang-Mills 
potential is basically abelian so that one has for the Wilson loop in any representation J: 

1 J i 

Wf = Wj M = —— £ exp im$, $ = I dxM? = / dS, - y (93) 

m=—J ^ 

In the last equation we have used the normal Stokes theorem for the circulation and also 
the fact that in the string gauge the magnetic field of the monopole is the Coulomb field of 
a point charge; dSi is the element of the surface spanned on the contour, and is orthogonal 
to the surface. 

Eq. (p3| ) is the gravitational Wilson loop for arbitrary contours but in a specific metric 
given by eq. (|88|). In order to generalize it to a general metric flS9]) all one needs is to perform 
the general coordinate transformation of eq. fl93|). To this end it is convenient to use, instead 



of dSi, its dual dS^ such that dSi = dS^ k . Taking into account that under the general 
coordinate transformation x % —>■ y l {x) the contravariant vector transforms as V 1 —>■ V k dky l , 
and the antisymmetric contravariant tensor transforms as dS^ — > dS mn d m y l d n y\ we get 
for the flux in eq. ( j93|) : 

$ = [dSi ^ = / dS** e ijk ^ - / dS mn e v kd ™f^y° y\ (94) 
J r J r 6 J \y\ 6 

This equation takes a more symmetric form if one uses the external coordinates 

3 3 

$ = (I) 2 2 / dSk €abc ^ ^ ^ W °' ^ = R- (95) 

Eq. (|93| ) together with eq. ( P5|) is our final result for the gravitational Wilson loop in the 
cylinder S 2 x R space of constant curvature R. Through eq. (89) the Wilson loop implicitly 
depends on the metric. Let us make a few comments. 
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The parallel transporter should depend on the metric along the contour but not on 
the surface spanned on the contour as it can be drawn arbitrarily. Despite the surface 
form of the result it is indeed so due to the fact that 



d k (e abc e ijk d t w a d jW b w c ) =0. (96) 
Therefore, the flux in eq. fl95|) can be presented as a circulation of a certain vector. 



The flux (|)5D coincides in form with a well-known expression for the winding number 
of the mapping S 2 i— > S 2 . For a closed or infinite surface the winding number is 
normalized as 



1 /2\§ 



In \R 



J dS k e abc e ijk d iW a d jW b w c = Q = integer. (97) 



For small contours eqs. ( |9~3| , |9~5D reproduce the result of the previous section. To check 
it, let us rewrite the general small- loop expansion (|69|) for the concrete metric fl37l). 
We find: 

2 R 

Rklpq = ekiu x u e pqv x v , g ij = -r 2 5 ij . (98) 



Putting it into eq. ( |69[ ) and then performing a general coordinate transformation x l — > 
y l (x) we obtain after simple algebra 

W o ! AJ+}) ( ejiurMMA^ ) 2 (99) 

J 6 V \y\ 3 J 



which coincides exactly with the expansion of eq. (|93| ) in small loop areas A5* up to 
the second order. 



8 Non-Abelian Stokes theorem in d = 3 gravity 

In section 4.1 we have shown that the gravitational Wilson loop in integer representations 
J as a functional of the metric is equal to the Yang-Mills Wilson loop as a functional of 
the Yang-Mills potential, provided one takes this potential equal to the spin connection 
corresponding to the metric in question, see eq. (|37D. For holonomies with half- integer 
J only the representation via spin connection is available. In this section we present the 
holonomy both for integer and half-integer J in a unified way, by introducing the 'non- 
Abelian Stokes theorem' similar to that for the Yang-Mills Wilson loop. We shall show 
that, being written down in a surface form, holonomies for both integer and half-integer J 
are actually expressible through the metric tensor and its derivatives only: the vielbein and 
spin connection (not uniquely defined by the metric) are, in fact, unnecessary. 
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For any representation J we can use our non-Abelian Stokes formula ( fH|) where the 
Yang-Mills potential is directly related to the spin connection: 



Wy = Wj M [spin connection] (100) 
= J Dn5(n 2 - 1) exp y J d 2 S l] [-F^n a + e abc n a {D.jif (D jn ) c . 

We shall replace the element of the surface by its dual, dS 1 ^ = e % ^ v dS p . Our goal will be to 
rewrite this representation of the Wilson loop in terms of the metric of the curved 3d space. 
To this end we first decompose the integration unit vector n in the dreibein: 

n a = m' e", n a n a = m l m j e^e^ = m l m j gij = 1. (101) 

The new 3- vector m is a covariantly unit vector. Since the background metric is fixed we 
just change the integration variables from n to m, the new integration measure being 

j Dn5(n 2 - 1) ... = J Dm^5(m i m j gij - 1) . . . (102) 



We next use the relation (f42|) of the field strength Fh computed from the spin connection 
A°: = (l/2)e abc uj bc , to the Riemann tensor. The first term in the exponent of eq. (|101|) 



becomes 



Using 



V9 

eq. fllQ3Q can be continued as 



first term = -dS p e ijp (-^ e abc m n e a n R l kij e b e ck . (103) 



e abc e bl e ck = J_ e lkm ^ ^ = ^ g a ( 1Q4 ) 



fll ,tten„^V^^ £ - 9m „ m «. (105) 
The combination of the covariant Riemann tensor and two antisymmetric epsilons has been 



encountered in section 5: in d = 3 it gives the Einstein tensor, see eq. ([73]) . We get therefore 



first term = dS p y/g {R5 p n - 2R p n ) m n , (106) 

where R v n is the Ricci tensor and R = R^ is the scalar curvature. 

We now turn to the second term in the exponent of eq. (|101|) and again use the decom- 
position ( |101j ). We exploit the fundamental relation ( |33| ) which can be presented as 

Dfn b ' = e b k {V j ) k l m l (107) 

where D bb ' = dj5 bb ' + e bcb ' Aj is the Yang-Mills and (Vj)f = dj5f + is the gravitational 
covariant derivatives. Therefore, the second term is 

second term = dS p e abc e a k e b e c n e ijp m k (W i ) l l ,m 1 ' (Vj)^m n ' 

= dS p ^ge i ii>e klnm k (V i rn) l (V :j m) n . (108) 
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Gathering eqs. ( |102| , |106| , |108| ) together we get finally a non-Abelian Stokes theorem for 
the gravitational Wilson loop or the trace of the parallel transporter for spin J along a closed 
contour: 

Wf = J Dm y/g 5{m i m j g i} - 1) 

x exp i J - J dS k yfg [(r5* - 2R k ) m p + e l]k e pgr m p (V, mf( V,- m) r ] . (109) 



Several comments are in order here. 

• The holonomy being defined as a path-ordered exponent is expressed here by a simple 
exponent of a integral over the surface spanned on the closed contour. That is why we 
call it a 'Stokes theorem'. The price to pay is a functional integration over covariantly 
unit vector m defined on the surface. 

• Eq. ( |109| ) is invariant under diffeomorphisms, in the sense that if one makes a general 
coordinate transformation x l —>■ x n (x l ) and changes the surface appropriately, the 
holonomy remains invariant. 

• The parallel transporter depends only on the contour but should not depend on the 
way one spans a surface on that contour. The surface integral in eq. ( |109| ) has the form 

JdS k ^gV k , (110) 
and the condition that it does not depend on the form of the surface is 



dk(VgV k )=o, (in) 

being equivalent to the condition 

(y k )iv l = o, (ii2) 

since T k t = Y k k = di In yfg. The check of eq. ( (L12j) is rather lengthy and we relegate it 
to the Appendix. 

Condition ( |112j ) or equivalently ( |111| ) being satisfied means that the surface integral 
can be written as 

J dS k yTg V k = J dS k e^%B k = -jdx l B, L (113) 

proving that it depends only on the contour, as it should be. However, the vector field 
Bi cannot be uniquely determined from the metric tensor and the covariantly unit 
vector m. 

The following comment is closely related to the previous. Parallel transporters of 
integer spins 1,2,... are defined via Christoffel's V symbols and hence by the metric 
tensor, while parallel transporters of half-integer spins 1/2, 3/2,... are not: they are 



19 



defined via the spin connection which is not uniquely constructed from the metric. 
Nevertheless, it should be expected that the holonomy for half-integer spins, i.e. the 
trace of a parallel transporter along a closed loop, being a diffeomorphism-invariant 
quantity, should be expressible through the metric only. The above eq. ( |109| ) solves 
this non-trivial problem: only the metric and its derivatives are involved. 

The solution is possible only when one presents the holonomy in the form of a surface 
integral, as in eq. (|109|) . One cannot do it in a contour form as it is not uniquely 
expressible through the metric. Were that possible, one would be able to write down 
a parallel transporter in terms of metric along an open contour as well, but that is not 
so for half-integer spins. 

Eq. ( |109| ) solves another long-standing problem, this time in the Yang-Mills theory. 
In another paper [TT|] we have shown that the SU (2) Yang-Mills partition function in 



three dimensions can be exactly rewritten in terms of gauge-invariant quantities which 
happen to be the six components of the metric tensor of the dual space. The usual 
argument why such rewriting is not too useful is that external sources couple to the 
Yang-Mills potential and not to gauge-invariant quantities. However, now we have 
demonstrated that a typical source, i.e. the Yang-Mills Wilson loop can be expressed 
not only through the potential but also through the metric tensor which is gauge- 
invariant. Thus, not only the partition function but also Wilson loops in the d = 3 
Yang-Mills theory can be expressed through local gauge-invariant quantities. 



9 Non-Abelian Stokes theorem in d = 4 gravity 

The aim of this section is to present the holonomy j \ in curved d = 4 space in the 
representation (Ji, J2) in terms of the metric tensor and its derivatives. Eq. fl6"8|) presents the 
holonomy in terms of the (anti) self-dual parts of the spin connection not uniquely determined 
by the metric, which is not satisfactory. In addition, we would like to get rid of the path- 
ordering in the Yang-Mills Wilson loops W n,p entering eq. fl68l). Both goals are achieved via 
the non-Abelian Stokes theorem similar to that of the previous section, which we are going 
to derive now. 

We start by applying the representation (|TJ]) to the Yang-Mills Wilson loop W": 

WJ = J Dn5(n 2 - 1) exp d J dS» v [-F a ^)n a + t ahc n a {D^)n) b (A,(vr)n) c ] (114) 

where D ab (-K) = 5 ab + e acb is the covariant derivative with respect to the self-dual part 
of the spin connection and F^ u (tt) is the appropriate field strength (|53"D ; it is related to the 



Riemann tensor via eq. (^). Let us introduce an antisymmetric tensor 

m KA = - n a T] aAB e AK e BX . (115) 

The first term in eq. ( |114j ) can be written as —R K \^ u m KX . The tensor m KA has actually only 
two independent components. To see this we introduce two covariant projector operators, 

= \v aAB V aCD e A e B e c ^ = ^(g K ,g Xu ~g KU g x , + ^e KXllu ), (116) 
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P K x^ = ^V aAB fj^ D eiefe^ = ^{g^g Xv -g Kv g^-^ge KXfMU ), (117; 



satisfying projector conditions, 

p± w' r>± _ p± (na} 

P^9^9 VV 'P^ pa = 0, (119) 
P&riTij* = 3. (120) 

-freAiii/ are (covariantly) orthogonal projectors, each having 3 zero and 3 nonzero eigenvalues. 
They project a general antisymmetric tensor into (covariantly) self-dual and anti-self-dual 
parts, respectively. It is easy to check that the tensor m nX introduced in eq. (|115|) is self-dual, 



and satisfies the normalization, 



P^m HX = 0, (121) 



m KX m KX = P+ ^ m KA wT = 1, (122) 

being a consequence of the normalization n 2 = 1. Therefore, m KX has, indeed, only two 
independent degrees of freedom in a given metric. We change the integration variables in 



eq. ( |114| ) from n to 



ni 



J Dn5(n 2 -!)••• = / Dm* x ^g S^P-^-mT) 5(m KX m KX - 1) . . . (123) 
Let us now compute the covariant derivative of m KX : 



= d^m KX + Tl v m uX + T x u m KU 



1 „aAB 



2* 



,a „Ak „B\ i „a/Q „Ak i p/t „Av\„BX i „a„i(t/o „BX , -pA „Bv^ 

\ aAB 
2 ' 



9/e i AK e BX + n a (<9 M e AK + T K ^e A,/ )e BX + n a e AK (d^e BX + Y x y 
d„n a e AK e BX - n a u AC e° K e BX - n a e AK u BC e cx 



(124) 



where in the last equation we have used the fundamental relation (|33|). We now insert the 
decomposition of the spin connection u AB into self-dual and anti-self-dual parts, eq. (f)T[). 
Using the relations for the 7], f] symbols, 

^AB^bAC = 6 ab 6 BC + e a b c r] cBC^ -aAB -bAC = gab gBC + ^c-cBC ? (125) 

v aAB-bAC = v aAC-bAB^ (126) 

it is easy to see that only the self-dual piece of ui AB survives in eq. ( |124j ) , giving 

m ^ = \rf AB e A * e BX (fy 5 ab + e acb irfy n b = hq aAB e AK e BX (D |4 (tt) nf . (127) 

In other words, the gravitational covariant derivative of m KX is expressed through the Yang- 
Mills covariant derivative of the n field, which is encountered in the second term of eq. (|114|) . 
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Using consecutively eqs. ( |125| , |127j) we obtain after simple algebra the final expression for 
WJ ( |114|) but this time presented in terms of the metric: 

x exp J dS^ [-i^ m KX - I ^ge KpaT g\\> m\ m%l . (128) 
Similarly, W p is obtained by integrating over anti-self-dual covariantly unit tensors: 

.J 2 



x exp *y J dS^ -R^ u m KX + -^g t Kpar g X \> m KX ' m\ m aT . v 



(129) 



As derived in section 4.2, the gravitational holonomy in representation ( J 1; J 2 ) is the product 
of the two, 



W 5uJ 2 ) = W l W Jf ( 13 °) 

Eqs. (|128| , |129| ) and (|130|) form the 'non-Abelian Stokes theorem' for the holonomy in curved 
d = 4 space. It expresses the holonomy via surface integrals spanned on the contour, and 
presents it in terms of the metric tensor and its derivatives only, without reference to the 
spin connection, even for half-integer representations Ji, J 2 . 



10 Conclusions 

The main result of this paper are the non-Abelian Stokes theorems for the holonomies: the 
Yang-Mills Wilson loop (eq. ( fill) ) and the traces of parallel transporters in curved d = 3 
(eq. (|109|) ) and d = 4 (eqs. (|128| , |129|) ) spaces. In all cases the path-ordered exponents of 
the connections are removed and replaced by ordinary exponents of surface integrals which, 
however, do not actually depend on the way the surface is spanned on the contour. The 
price to pay for the removal of path ordering is high: we obtain functional integrals instead. 
In the simplest case of the SU(2) Yang-Mills theory it is an integral over a unit 3-vector n 
'living' on the surface, in the case of the d = 3 Riemannian manifold it is an integral over a 
covariantly unit 3-vector m, in case of d — 4 one integrates over (anti) self-dual covariantly 
unit tensors. 

In spite of functional integration we believe that our formulae are aesthetically appealing. 
As compared to path-ordered exponents they are better suited to averaging over quantum 
ensembles of Yang-Mills fields or over various metrics. We hope that elegant formulae can 
be also used in more general settings. 

In addition to the general non-Abelian Stokes formulae we have presented holonomy as 
a surface integral for a specific background, namely for a d = 3 space of constant curvature 
with cylinder topology S 2 x R. The 'gravitational Wilson loop' is given by a formula for the 
character whose argument is the winding number of external coordinates, see section 8. 

Parallel transporters of integer spins have a dual description: one can define them either 
as a path-ordered exponent of Christoffel symbols or as a path-ordered exponent of spin 
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connection in appropriate representation. In section 4 we have shown that both representa- 
tions are equivalent. Even though spin connection is not uniquely determined by the metric 
tensor, this equivalence means that the holonomy written in terms of spin connection can 
be in fact expressed through the metric only. 

For half-integer spins the situation is far less trivial since the only way to define the 
holonomy is via the spin connection, and it is not at all clear beforehand that it can be 
uniquely written through the metric tensor and its derivatives. The non-Abelian Stokes 
theorem of this paper demonstrates that it is indeed so, but only when one presents the 
holonomy in the surface form which is uniquely determined by the metric. Though the 
surface integral does not depend on the way one draws the surface and can be actually 
written as an integral along the contour, the contour form is not uniquely defined by the 
surface one, and it reflects the ambiguity in determining the spin connection from the metric. 

This finding has interesting implication for Yang-Mills theory in three dimensions, which 
can be identically reformulated as a quantum gravity theory, with the partition function writ- 



ten as a functional integral over the metric tensor of the dual space ||11|| . This metric tensor 
is local and gauge invariant (in the Yang-Mills sense). However, one might wish to calculate 
the average of the Wilson loop which, originally, is defined by the Yang-Mills potential, but 
not by the metric tensor. In the 'quantum gravity' formulation the Yang-Mills Wilson loop 
becomes a parallel transporter in the gravitational sense. Therefore, it is very important 
that it, too, is expressible through the gauge-invariant metric tensor, in any representation. 
Thus, not only the partition function but also the Wilson loop can be presented in terms of 
local and gauge-invariant quantities. This subject will be presented in more details elsewhere. 

One of us (V.P.) thanks NORDITA for kind hospitality and the Russian Foundation for 
Basic Research for partial support, grant RFBR-00-15-96610. 

Appendix. Proof that eq. ( |tO§l ) does not depend on the 
surface 



The path-integral representation for the 'gravitational Wilson loop' ( |109|) should not depend 



on the choice of the surface but only on the contour on which the surface is spanned. To 
prove that it is indeed so, one has to check eq. (|112|) , 



(v*)fy' = o, (131) 

where 

V k = (R5 k p - 2R k ) mP + e ijk e pqr m p (V; m) q {V 3 m)\ m i m i g ij = 1. (132) 



To simplify notations we denote covariant derivatives by ";" (see |16||). Explicitly, the 
covariant derivatives of a scalar, vector, tensor are given by 

S-k = dk S, 

v; k = d k v* + r kl v\ v i;k = d k v i -v l ik v h 

T-i = d k Ti + V kl T l i + r{ l T i \ T.j.k^dkT.j-rl^-T'Tu, etc. (133) 
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An ordinary derivative of a convolution of two tensors can be written as a sum of covariant 
derivatives: 

<h, (v;. 1 ' r:r) = r^ k A t { }- + v;; l: -' T^: k . (134) 

We apply the covariant derivative to the first term of the vector V k : 

= (R5 k p - 2R k ) ± m p + (R5 k p - 2R k p ) m%. (135) 

The covariant derivative of the Einstein tensor is known to be zero ([0, eq.(92.10)). There- 
fore, only the second term in eq. ( |135|) survives. 

Next, we apply the covariant derivative to the second term of the vector V k : 



V, 



\R5 k — 2R k 



vtv 



Ajk 



~pqr 



m ? '(V i m) 9 (y 7 -m) r 



e ijh e 



pqr 



(V k m) p (y i m) q {Vj m) r + 2 e ijk e pqr m p (V 4 m) <? (V fc V i m) 



(136) 



The first term here is zero, for the following reasons. We differentiate the condition that m' 
is a covariantly unit vector, 







d k (m l m?gi^ = 2g^ (V fc m) 1 m? = 2 (V fc m)* m„ 



(137) 



since the covariant derivative of the metric tensor is zero. It means that three vectors, 
(^1,2,3 m T are n °t linearly independent as there cannot be three linearly independent vectors 
orthogonal to some vector (in this case m;) in three dimensions. The first term in eq. ( |136| ) 
is an antisymmetrized product of these three linearly-dependent vectors, therefore, it is zero. 

As to the second term in eq. ( |136j ) we notice that it contains the commutator of covariant 
derivatives, equal to 



1 



1 



e^CVfcVj m) r = -e^ k [V^l m s = -e^ k g rt R tskj m s (138) 

where Rtskj is the Riemann tensor. Therefore, the second (and the only nonzero) term in 
eq. ( |136j ) can be written as 



e ijk e pg r 9 rt Rtskj m p m s (V< m) 



(139) 



We next use eq. ( |70[) to express the Riemann tensor through the Ricci and metric tensors 
and write down the product of two epsilons as a determinant made of Kronecker deltas. 
Performing all convolutions we obtain that eq. (|139| ) can be identically rewritten as 



9qs 



2R 



2R 



m p m s ( Vj m) 



■g ps (R5' l q —q 

Here the first term is zero because of eq. (|137|) while in the second term we use g ps m p m' 
As a result we get 



2R\ 



(^m) 



(140) 
1 = 1. 

(141) 



which cancells exactly with eq. ( |135| ). Thus, (V&)* V 1 = 0, q.e.d. 
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